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Abstract 

Within the LAGUNA design study, seven candidate sites are being assessed for their fea- 
sibihty to host a next-generation, very large neutrino observatory. Such a detector will be 
expected to feature within a future European accelerator neutrino programme (Superbeam 
or Beta Beam), and hence the distance from CERN is of critical importance. In this article, 
the focus is a ^^Ne and ^He Beta Beam sourced at CERN and directed towards a 50 kton 
Liquid Argon detector located at the LAGUNA sites: Slanic (L = 1570 km) and Pyhasalmi 
(L = 2300 km). To improve sensitivity to the neutrino mass ordering, these baselines are 
then combined with a concurrent run with the same flux directed towards a large Water 
Cerenkov detector located at Canfranc {L = 650 km). This degeneracy breaking combina- 
tion is shown to provide comparable physics reach to the conservative Magic Baseline Beta 
Beam proposals. For ^^Ne ions boosted to 7 = 570 and ^He ions boosted to 7 = 350, the 
correct mass ordering can be determined at Slanic for all 6 when sin^ 26'i3 > 4 ■ 10"'^ in this 
combination. 
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1. Introduction 

Results from a series of atmospheric [H, solar jsj-Isj, reactor and long baseline 
accelerator jol, neutrino experiments indicate that neutrinos are both massive and mix 
amongst themselves. A combined analysis points to two approximate 2-neutrino mixing 
schemes, each parameterised by a mixing angle and a mass-squared splitting. The extent 
to which these two schemes combine into a single 3-neutrino picture is controlled by the 
size of a third mixing angle, 613. Defining Am|^ 
available data returns the best fit values 



rrij — mf, a combined analysis of all 



|Am3i| 
Amj^ -- 



= 2.47 X 10"^ eV^ and 
7.59 X 10"^ eV^ and 



sin^ e. 



23 



sin^ 9i2 



0.463 
0.319 



The third mixing angle, ^13, is constrained to be 

sin^ ^13 < 0.016 (0.053) at la {3a); (1) 
although some collaborations report a hint at low significance from a combined analysis of 



atmospheric, solar and long-baseline reactor neutrino data [12|, |13 



The long term targets for the long baseline neutrino oscillation programme is to determine 
the unknown neutrino mixing angle ^13; whether there is CP- violation in the lepton sector, 
that manifests itself through a phase 6; and to determine the sign of the atmospheric mass- 
squared splitting. The extraction of both the CP-phase and sign of the atmospheric mass- 
squared splitting require a three- neutrino analysis, the extent of the effect controlled by 6*13. 
Therefore sensitivity, optimisation and general experimental strategy are intrinsically related 
to the size of ^13. Running and near future experiments will be the first to probe ^13 below 
the current experimental limit. If ^13 is beyond the reach of these experiments, then intense 



Neutrino Factories 17, 18 



sources of neutrinos from next generation Superbeams |14|-|16 
or Beta Beams 19|] will be necessary. 

A Beta Beam is a variation on the Neutrino Factory that instead sources neutrinos from 
the decays of boosted radioactive ions. Proposed by Zucchelli in 2001 19|, a Beta Beam 



distinguishes itself by sourcing beams of z/g or z/g without contamination from other flavours 
and CP-conjugates. Through the analysis of the z/^ and z/^ appearance channels, the Beta 
Beam provides a competitive physics reach on the unknown elements of neutrino mixing 
that have yet to be extracted from experiment. A Beta Beam will be a facility that will 
exploit existing, or proposed, ion production facilities and accelerator infrastructure. The 
major addition required is a large decay ring to store the ions once they are boosted to 
the desired energy. In the laboratory frame, the neutrino flux is a function of both the ion 
decay Q-value and the Lorentz boost of the ions. Low Q-value ions; ^^Ne and ^He; and high 



Q- value ions ^B and ^Li have been identified as the best candidate ions 20|. The focus of 
phenomenological study has been on a Beta Beam source at CERN using either the Super 
Proton Synchrotron or an enhanced 1 TeV version, as required in some LHC upgrade 

scenarios |26m33i] (Other facilities could be the Tevatron at Fermilab 34j-l36| . or a re-fitted 
HERA- ring |37|.) In addition to specific setups, a number of optimisation studies have also 
been carried out 138-40 . 
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In this article, the physics reach of Beta Beams directed towards Liquid Argon detectors 
in Europe will be simulated on the assumption that a 1 TeV machine will be available (4ll |. 
The study of both low and high boost Beta Beams directed at large Water Cerenkov (WC) 



detectors has been studied in detail and is well understood j21|-|27l|. Water Cerenkov de- 
tectors are best suited to short basehnes since they only use the quasi-elastic events to 
reconstruct the signal. For longer baselines, and hence higher energies, it is necessary to use 
and reconstruct the energy of multi-particle final state interactions, especially deep inelastic 
scattering events. For this reason, one must use active scintillator, calorimeters and projec- 
tion chamber technologies for intermediate and very long baselines. This study considers 
the GLACIER liquid argon detector j42| as a far detector for the baselines CERN-Slanic 
(L = 1570 km) and CERN-Pyhasalmi (L = 2300 km) for 7 = 350, 350 and 7 = 570, 350 ^^Ne 
and ^He Beta Beams. (The CERN-Boulby baseline with a non-WC detector was discussed 
in jlil for different exposures and so the very similar baseline for Sierozsowice {L = 950 km) 
will not be included here.) These two baselines will then be combined with the same neutrino 
flux directed towards a large Water Cerenkov detector based at Canfranc {L = 650 km) to 
examine the possible broad physics reach. 

This work is part of the LACUNA design study: an EC-funded project to assess the 
feasibility of underground laboratories in Europe capable of housing a future, large neutrino 
observatory. In addition to long baseline physics, the detectors aim to improve the bounds 
on proton decay half-lifes up to ~ 10^^ years to test a number of theoretical models; to 
detect neutrinos from astrophysical objects and Cosmological sources; continued study of 
solar and atmospheric neutrinos; and new sources such as from Dark Matter annihilations in 
the Earth's centre or Sun's core. Work towards understanding detector response to various 



neutrino sources in currently on-going. This work is one of a series of studies |16l . l29l . 143 
examining the physics reach on long baselines experiments with the LACUNA network of 
detectors based on available information. These studies will help prioristise the potential 
host laboratories for the future neutrino observatories. 

The remainder of this article is organised as follows; in Sec. El the baselines available 
within Europe as part of the LACUNA design study [441] will be listed and the basic phe- 
nomenological strategy for short and intermediate baselines will be summarised. In Sec. |3] 
a description of the simulations carried out is made followed by the results for Beta Beams 
directed along single baselines. In Sec. lU the simulations from Sec. E] are augmented with 
an additional Beta Beam directed towards a large Water Cerenkov based at Canfranc. The 
results are then discussed in relation to Magic Baseline Beta Beams in Sec. [51 Finally, in 
Sec. [6l the study is summarised. 

2. Laguna sites hosting Beta Beam far detectors 



Detector R&D is not part of the LACUNA Design Study [4^ ; however the physics reach 
of these sites with respect to any European accelerator neutrino programmes will be critical 
in any decision making process. Since future European Superbeams and Beta Beams will 
be one of the main users for these detector technologies; the distance from CERN will be 
a critical factor in the final decision. However, since the size of 613 and the number of 
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future facilities and their scale is unknown, it is very hard to determine what is the optimal 
baseline and detector choice. More precisely, optimisation of a facility can only be carried 
out once its purpose and circumstance is defined. Since ^13 is unknown, and (specifically) 
the role in which a Beta Beam will take in a long term experimental strategy is uncertain; 
optimisation of the Beta Beam is not a well defined process. A typical strategy, therefore 
is to aim for a broad physics reach using the ability to rule out 6*13 = 0°, 6 = 0° and 180°, 
and determine the correct neutrino mass ordering as indicators for a particular experimental 
setup; an approach also adopted here. How one achieves this is dependent, amongst other 
things, on the choice of baseline (or baselines) and detectors technologies. 

Within Europe there are three detector options being considered: a single volume Liquid 
Argon time projection chamber (LAr) known as GLACIER [42]; a 500 kton Water Cerenkov 
(WC) known as MEMPHYS |45| suitable only for energy reconstruction at short baselines 
where quasi-elastic events dominate; and a non-segmented liquid scintillator detector known 



as LENA |46[. The 7 laboratory sites being considered within LAGUNA to house these 
detectors are listed in Tab. [T] along with their distances from CERN and corresponding first 
oscillation maximum energies. It clearly seen that a wide range of neutrino baselines are 
possible in Europe ranging from L = 130 km to L = 2300 km. 

To develop a strategy, or interpret results of numerical simulations, for particular base- 
lines it is useful to consider an analytical approximation to the oscillation probability. With 
the definitions a = and A = -^^-^ ; the assumption of a constant matter density 

profile along the baseline A = ^^^a""'^ ; and by expanding in all the small parameters one 



finds 47 



where 



P{u, ^ u^) ~ sin^ 2^13 ■ Ti + a ■ sin ^13 ■ (T2 + T3) + ■ , (2) 
sin2[(l - A) ■ A] 



Ti = sin^ 6123 

sin(AA) sin((l - A) A) 



(l-Ay 
■ sin 6c p ■ sin 26'23 ■ sin A' 



T3 = cos 6c p ■ sin 26'23 " cos A 



A (l-A) 
sin(AA) sin((l-A)A) 
(l-A) 



T4 =cos 6'23 ■ sm 26'i2 ; (3) 

As is well known, and clear from this expression, the determination of and 6 is severely 



affected by parameter correlations and degeneracies [48|. For a given bin at fixed baseline, 
up to 8 possible parameter sets can fit the data. The challenge for a future long baseline 
experiment is to successfully resolve these degeneracies and push for a good physics return 
over the sought region of parameter space. With a 1 TeV machine, it is in principle possible 
to source a Beta Beam from very short long-baselines (L = 130 km) up to very long long- 
baselines (Magic Baselines j33[). A Beta Beam proposal therefore has access to two types 
of baseline in which some of the degeneracy can be naturally suppressed: 
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Location 


Distance from 


Energy 1st Osc Max. 




CERN [km] 


[GeV] 



Frejus (France) 


130 


0.26 


Canfranc (Spain) 


630 


1.27 


Umbria (Italy) 


665 


1.34 


Sierozsowice (Poland) 


950 


1.92 


Boulby (UK) 


1050 


2.12 


Slanic (Romania) 


1570 


3.18 


Pyhasalmi (Finland) 


2300 


4.65 



Table 1: Potential sites being studied with the LAG UNA design study [4j|. The Umbria site is of interest 
if an upgrade to the CNGS beamline is pursued. 

1. At short long-baselines (L < 700 km say), the matter effect is sufficiently small so that 
the true {613, S) and the fake solution corresponding to the incorrect mass ordering 
are close together. Consequently, the sensitivity to 613 and 6 is typically very good, 
especially with the availability of large Water Cerenkov detectors. Using the above 
expression, the CP- violation contribution to the probability is maximal for L/E = 
515/GeV. A Beta Beam flux with Lorentz boost 7 = 350, matches well the CERN- 
Canfranc baseline {L = 650 km) [3o| . 

2. At the 'Magic Baseline', cancellations leave only the atmospheric contribution (Ti) to 
the appearance probability and thus there is no S dependence. Numerically, this is 



found to be at L = 7250 km [49||. Use of the Magic Baseline in isolation will therefore 
return a reasonably clean measurement of 6*13 and sign(Am3]^). The excellent sensitivity 
to these parameters, which is typical to proposals incorporating it, is because of the 
proximity to the matter resonance which, with the larger cross-sections at high energy, 
can compensate for the dependence of the un-oscillated neutrino flux. 

Both these options are a staple in long baseline proposals since the effective removal of 
one of the unknown parameters helps greatly with degeneracy resolution. However, a extra 
baseline will be always be needed if one is to have a competitive reach on all three physics 
indicators: discovery of non-zero ^13, CP non-conservation, and ability to rule out the in- 
correct neutrino mass ordering. For Beta Beams, the optimal choice for a second detector 
has not be determined rigorously, although several possibilities have been put forward for 



accompanyment of the Magic Baseline |30l-l32 



In this article, a Liquid Argon detector is considered at the two longest baselines within 
LACUNA and therefore there is no natural suppression of the degeneracy. By itself, a single 
baseline Beta Beam using a 50 kton detector will not return competitive sensitivities at the 
longer LACUNA baselines owing to a suppression of the neutrino flux from the baseline 



(oc -L~^) and the presence of parameter degeneracies 48|. In isolation, binning the data to 
extract the oscillatory structure of the appearance probability is the required strategy to 
break parameter degeneracies. The relative weight of the atmospheric (Ti), interference (T2 
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and T3) and solar (T4) contributions to the probability changes with the neutrino energy for a 
fixed baseline. Since different contributory features of the appearance probability dominate 
different bins and that the location of parameter degeneracies are energy dependent, this 
strategy greatly aids degeneracy resolution. 

The above strategy was discussed in [28i] for the CERN-Boulby baseline {L = 1050 km) 
and is used in a number of other long baseline proposals, such as the Wide Band Super- 
beam [50!] and the low energy Neutrino Factory 18 . It was argued semi-analytically that 
combining data from first and second oscillation maximum, even for a neutrinos (or anti- 
neutrinos) only, helps considerably in breaking the intrinsic energy degeneracy (or energy 
degeneracy [5l|). Degeneracy in S, but not ^13 still remains, however. Incorporating data 
from surrounding bins typically allows determination of the true 6. Resolving the sign(Am3]^) 
degeneracy at short to intermediate baselines is far harder. The difference in the probability 
for normal and inverted ordering results from two effects: the sin A in T2 (which is present 
even in the vacuum for 3-neutrino mixing); and from {1 — A) in Ti, T2 and T3. The discrep- 
ancy between normal and inverted orderings clearly increases with baseline and neutrino 
energy. For the baselines considered here, the approach is therefore to extract the sign from 
high energy bins. Since these bins also contain information on 6*13 and 6, degeneracy will 
need to be resolved. This, again, is achieved through the combination with the low energy 
bins where the probability splitting between mass orderings is much smaller and hence the 
degenerate solution is much closer to the true solution. The limit on the sensitivity to the 
mass ordering is then limited by the available neutrino flux. 

In the following section, the effectiveness of the above strategy will be explored for a 
Beta Beam directed along the CERN-Slanic and CERN-Pyhasalmi baselines. The event 
rates from the low energy bins are an important component of this approach, however, they 
are small owing the small detector mass and cross-sections. Therefore, following the initial 
analysis, the Liquid Argon baselines will be put in combination with a large Water Cerenkov 
based at the Canfranc laboratory. As indicated above, this will provide a clean measurement 
of 6*13 and 6 which can then be combined with the high energy bins of the Liquid Argon 
detectors. The reach on the correct mass ordering is expected to improve signiflcantly as the 
inclusion of the extra (short) baseline is akin to enhancing the event rate in the low energy 
bins. The analysis carried out will be discussed in more detail in the following sections. 

3. Single baseline study 

The physics strategy in this study exploits Ue and z/g beams sourced from the decays 
of boosted ^^Ne and ^He. With production and acceleration of the ions based at CERN, 
the physics will be simulated for a 50 kton Liquid Argon detector located at the Slanic 
{L = 1570 km) and Pyhasalmi {L = 2300 km) mines. In the flrst instance only these 
baseline will be considered for a 5 year run of Uf, and 5 years of z/g for the two boost 
pairings (7^,7^) = (350,350) and (570,350). These boost assignments make the standard 
assumption that a 1 TeV Super Proton Synchrotron will be available for a Beta Beam based 



at CERN |41[. In the following section, the simulations will also include the combination 
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with a larger Water Cerenkov detector based at Canfranc {L = 650 km) exposed to the 
same neutrino source. 

The current R&D in Liquid Argon detector technology is working towards a target mass 



of 100 kton [42|. However, in this article a 50 kton detector is used for the following reasons 



To bring it in line with other non-WC based Beta Beam studies in the literature. 
The sought level of ^^Ne decays along the straight section of the storage ring for the 



techniques [52 



7 = 100, 100 proposal |22| currently is beyond reach by a factor 20, using ISOLDE 



A 100 kton detector is assumed but with a factor of 2 used to offset 
a smaller than sought decay rate. 

The ion decay rates in the straight sections of the decay ring are currently unknown. R&D 
towards production of ^^Ne and ^He has focussed on the use of ISOLDE techniques within 
the EURISOL design study |53| with target rates of 1.1 x 10^^ and 2.9 x 10^*^ yearly decays 



for a 7 = 100, 100 machine for ^^Ne and ^He respectively. Currently, the projected rates 



are a factor 20 short for ^^Ne and about 2 short for ^He 52|. This in part is because ISOL 
techniques, as applied by Nuclear Physicists, are optimized for more exotic nuclides that are 
of little interest for Beta Beams [54] • The route to the sought decay rate for ^^Ne could be 



through 'direct production', for example through the ^^0(^He,n)^^Ne reaction 55 . 

Atmospheric neutrino events are skewed below 1 GeV and cause a problem for sourcing 
decays for the short baselines and need to suppressed with a restrictive duty factor in the 
decay ring (~ 10~^). However this constraint can be loosened for longer baselines where 
neutrino energies > 1 GeV are more important. This action will be necessary to reclaim 
the decay rate for high-7 machines which, to first order, scales as 7^^. (The impact of how 
much one inhibits the atmospheric neutrino background at the longer baselines has yet to 
studied in detail.) It may be possible to also claim some deficiency in the event rate through 
loosening the duty factor. To aid comparison with the latest studies, in this article the 
assumption that 3 x 10^^ useful ion decays per straight section of the decay ring can be 



sourced per year for each ion as suggested in |5J], and used in [30|, will be adopted. 



For short and intermediate baselines, the aim is to exploit the energy dependence of 
the oscillatory structure of the appearance signal to help break the parameter degeneracies 
and push for a good physics reach. This is achieved through binning the signal to separate 
out the low and high energy appearance events so that the different strengths of the solar, 
interference and atmospheric contributions to the appearance probability can be observed. 
To this end, the detector has been assumed to possess a low energy threshold of 0.4 GeV and 
has binned up to 2.0 GeV in 0.2 GeV intervals. All bins thereafter are in 0.5 GeV intervals 
up to the maximum laboratory energy of the neutrinos. At present, detector response for 
incident Beta Beam fluxes have not been simulated for the GLACIER Liquid Argon detector. 
Work in this respect in ongoing, with response data soon to be available for Superbeams 



directed at this technology 56|]. In its absence, a fiat event efficiency of 80 % has been taken 
for all channels with signal errors set at 2.5 % and an energy resolution of cr{E,^) = 0.15eJ^. 



^ These assignments are conservative with studies for FLARE [57| indicating that below 1 GeV, resolution 
could be as good as 2 %, whilst above 2 GeV a{Ey) ~ Q.lQEy. 
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Figure 1: 3cr C.L. contours for discovery for non-zero 6*13 (left), CP-violation (centre), and sign(ATO§]^) 
determination (right). In each plot, the solid line corresponds to a 7 = 350,350 ^®Ne and ^He Beta Beam 
for a 50 kton Liquid Argon detector located at Slanic {L = 1570) (top line); and Pyhasalmi (L = 2300 km) 
(bottom line). The dashed lines correspond to a 7 = 570, 350 ^^Ne and ^He Beta Beam. 



All simulations in this study have been carried using the GLoBES long baseline simula- 
tion software ^] . The known oscillation parameters have been fixed to their current central 



values, taken from [11|, and are always marginalised over. The exception is 623 which is 
fixed to 45° so that the octant degeneracy is absent. The errors have been set to 1 % for the 
solar parameters and 5 % for the atmospheric parameters and matter potential. Negligible 
background from neutral current events is expected and has been set to 0.1 %. The 1 degree 
of freedom convention is adopted for all sensitivity plots in this paper. 

3.1. Results 

The physics reach of the two baselines considered in summarised is Fig. [T] The 3cr 
confidence level contours are shown for sensitivity to non-zero 613 (left), CP non-conservation 
(middle), and to resolving the sign(Am3]^) ambiguity (right). The analysis takes into account 
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the impact of intrinsic and the sign(Am3^) degeneracies. The top hne of Fig. [T] shows the 
results for the CERN-Slanic basehne, whilst the bottom display the outcome of the CERN- 
Pyhasalmi simulations. 

The best sensitivity to non-zero ^13 and CP-violation is found for the CERN-Slanic 
baseline. This is to be expected since, with the same source, the flux is larger for this 
baseline. The weaker matter effect means that the sign(Am3^) degeneracy interferes less 
with these measurments. For the 7 = 570, 350 boost pair, non-zero di^ can be seen down to 
sin^ 26 ~ 10~^, and, for both boost pairs, at all values of 5 for sin^ 26'i3 > 10~^. However, 
there is a marked difference between the two boost pairings for sensitivity to CP-violation. 
For the 7 = 350, 350 pair, the ability to rule out 5 = 0° or 180° is restricted to sin^ 26*13 > 
10~^, but increasing the boost of the ^®Ne ions returns a large region of parameter space for 
5 < 0° and centred on sin^ 26*13 = 10^^. There is little enhancement on the region for low 
boost paring. This is suggestive that degeneracy is a problem for the 7 = 350, 350 boost 
pairing; especially for sin^ 2^13 ~ 5 • 10~^ — 1 ■ 10~^ where there is a gap in CP-violation 
sensitivity. The lower event rates for the longer baselines mean that the data is insufficient 
to reduce the significance of some degenerate solutions. Although, the ability to rule out 
the incorrect mass ordering is poor relative to Neutrino Factories, for a Beta Beam it is not 
intrinsically bad. For the high boost run, the correct ordering can be indentified down to 
sin^ 2^13 = 2 • 10~^, with determination for all values of 5 for sin^ 2^13 > 10~^ in both cases. 
Although the increase in the ^^Ne boost improves the reach, it does not do so significantly. 
Increasing the boost makes data from higher enegies available without significantly altering 
the event rate and composition at lower energies. Since European baselines make use of 
these bins in combination, improving one without the improving other need not, and has 
not, dramatically improved the physics return. The low event rate at low energies is still 
insufficient to break the degeneracy for small values of ^13. 

The physics return for the CERN-Pyhasalmi baseline is weaker for each of the physics 
indicators, with little sensitivity for sin^ 26*13 < 10^^. Principally, this is due to the 
dependence of the un-oscillated neutrino flux. In particular, the ability to determine the 
correct mass ordering is very poor even given the large matter effect at this baseline. The 
true and incorrect mass ordering solutions will be sufficiently separated in (6*13, 5) space; 
however, the low event rate means that the solution regions at 3(j will be large and possibly 
merged together. When combined with a large solution region from the low energy bins 
(also owing to low event rates), the data is insufficient to break the degeneracy for small ^13. 

In summary, neither of the baselines considered here have, in isolation, the right com- 
bination for degeneracy breaking ability and sufficient un-oscillated event rate to return 
a competitive physics reach on all physics indicators. One solution to this problem is to 
combine these baselines with another beam. It has been proposed in 35|, |59|| to use combi- 
nations of Beta Beams and Superbeams at these shorter baselines to break degeneracy and, 
in particular, improve the ability to determine the mass ordering. This strategy exploits 
the different forms of the appearance probability for the i^ti, i^e ^ti,i^ti and 

channels which in turn affects the location of the degenerate solutions for different 
oscillation channel pairs. 

In the next section, a variant on this approach will be adopted. Specifically, a second 
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Beta Beam baseline will be introduced in view of returning sensitivity to ^^13 and 6 to much 
greater precision. This is extrapolated from the strategy advocated earlier of using low and 
high energy bins and is motivated from the results for the Slanic baseline indicating that 
improvement in the high energy neutrino event rate is insufficient to significantly improve 
the hierarchy reach. Increasing the ion boost does not significantly improve the event rate 
at low energies; however, introducing a second Beta Beam at a shorter baseline, and with 
a large Water Cerenkov will have the desired effect. This shorter baseline can be used to 
locate the true solution accurately to the level sin^ 29i3 ~ 10~^, effectively replacing the low 
energy bins of the (now) far detector. With this information the longer baseline can then 
rule out the incorrect solution from the sign(Am3^) degeneracy for much smaller sin^^is. 
This idea will be presented in the next section. 



4. In combination with a large Water Cerenkov 

As discussed in Sec. [21 for a single baseline, the strategy is to use the complementary 
information from low and high energies to break degeneracies and push to physics reach 
at small 613. However, for the intermediate baselines considered here, owing to smaller 
detector sizes and increasing sign(Am3]^)-degeneracy dominance, the ^13 and 6 reaches are 
relatively poor compared to the Water Cerenkov Beta Beams. This was demonstrated 
through the simulations of the previous section which indicated that degeneracy persisted 
for the intermediate baselines of Europe. A straightforward way to overcome this is to 
increase the exposure so that degenerate solutions lose statistical significance. However, for 
the Beta Beam, demands on production and ion storage impose restrictions on scaling up 
of this kind. 

The challenge is improve the event rate at the low energies and then use this data in 
combination with high energy neutrinos. This approach is distinct from the standard tech- 
nique to resolve the neutrino mass ordering. It is typical to consider very long baselines 
where the matter effect is strongest and the closeness to the matter resonance allows the 
recovery of flux reduced by its dependance. Most interest centres on the Magic Baseline 
where CP-violation effects vanish from the appearance probability, leaving only the atmo- 
spheric contributions. There have been a number of studies using a second Beta Beam in 
addition to a Magic Baseline Beta Beam to help break degeneracies. Thus far, a Magic 



Baseline/short baseline combination has been studied in 30|, |3l| and a Magic Baseline/long 



baseline combination in [32|. The strategy employed in this section is to deliberately place 
the far detector at a baseline where the degeneracy is strong. In a short baseline/intermediate 
baseline conflguration the short baseline will provide the sensitivity to ^13 and 6 down to the 
sin^ 2^13 ~ 10~^ level. This information can then used to rule out degenerate solution regions 
in (^13,5) space for the longer baseline. (For a short baseline, the sign(Am3^)-degenerate 
solution will be very close to the true solution. As the matter increases, the two solutions 
diverge.) The combination of data sets is expected to improve the sensitivity considerably. 
This is demonstrated in the following simulations. 
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Figure 2: 3a C.L. contours for discovery (for the far detector at Slanic) for non-zero 6*13 (top left) and 
CP-violation (top right). On the bottom line 3(j C.L. contours for sign(Am|]^) determination for Slanic 
(left) and Pyhasalmi (right) In each plot, the solid line corresponds to a 7 = 350,350 ^*Ne and ^He Beta 
Beam and the dashed lines correspond to a 7 = 570,350 ^^Ne and ^He Beta Beam. 



4.1. Results 

In Fig. 121 the 3a confidence levels for the three physics indicators are shown for a large 
Water Cerenkov, based at Canfranc, in combination with the same Liquid Argon detector 
as previous for the the Slanic and Pyhasalmi laboratories. The Water Cerenkov dominates 
the and 6 sensitivities and the reaches are essentially identical in both combinations. 
Therefore only the results for Slanic are shown here. The CP-discovery plot is both smooth 
and symmetric. There is no residue intrinsic degeneracy at sin^ 29i3 ~ 10~^ as its location is 
different for each baseline and therefore can be resolved. For both baseline pairs, non-zero 
^13 and S distinguishable from 0° and 180° can be achieved down to sin^ 2^13 = 5 • 10^^. 

The important result is that the combination of baselines does indeed improve the ability 
to rule out the incorrect mass ordering beyond the capability of either baseline separately. 
Specifically, for 7 = 350, 350 and for both combinations, the correct mass ordering can be 
extracted for all 6 down to sin^ 29i3 ~ 6 — 7 ■ 10^^. For 7 = 570, 350 this improves to 
sin^ 2^13 ~ 4 ■ 10"'^ which is the level reported for the conservative Magic Baseline proposal 
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in [30|. Even for 7 = 350,350, there is substantial resolving power at the 4 ■ 10"^ level 



indicating that a more more minimal Beta Beam is capable of returning similar physics. 
5. Discussion 

The simulations presented in this article indicate that a two-baseline Beta Beam con- 
figuration using a short and intermediate baseline can obtain a similar physics reach as 



conservative Magic Baseline proposals, such as in [30[. Long baseline experiments incorpo- 
rating the Magic Baseline aim to exploit the absence of the CP-violation at this baseline 
along with the nearby resonance associated with atmospheric neutrino mixing. Since such 
an environment provides a relatively clean measurement of the mass ordering, usually to 
very small ^13 , Beta Beam studies aiming to achieve the best reach in this respect typically 



exploit it |30H33|. It is instructive to highlight the reason why a more minimal setup can 
return similar sensitivity to the conservative Magic Baseline proposals. 

The crucial point to realise is that the sought features of the baseline, namely no CP- 
phase and proximity to the matter resonance, are features of the appearance probability. 
In an experiment, one measures the event rate which is a convolution of the unoscillated 
neutrino flux at detector, a cross- sect ion, detector efficiencies and the probability. Having 
a signal clean from CP-violation is of little use if the event rate is too low to provide a 
competitive sensitivity. The unoscillated flux has a dependence which greatly reduces 
the flux at very long baseline. This reduction can be recovered through the larger cross- 
sections at high energy and the matter enhancement of the appearance probability. On the 
assumption that the same number of ion decays is available for the very long baselines as 
for the short, the number of appearance events is roughly independent of the baseline at 



very long baselines [33[. In this instance, the ability to extract the mass ordering will peak 
where the signal is cleanest from 6 'contamination': the Magic Baseline. However, a second 
baseline is always necessary to have access to the CP-violation. In which case, it may be 
possible to use synergy between baselines to extract the sign of the mass splitting at a similar 
level as demonstrated in this study. 

Never-the-less, the sensitivity presented here can still be beaten by most Magic Baseline 
Beta Beams considered. The point raised in |30i] is that construction restrictions on the size 
of the decay ring could result in a considerable reduction of the flux at very long baseline 
in a realistic proposal. (Longer baselines require the decay ring to dip at larger angles to 
the surface so that the maximum depth of the ring could be beyond 2km if very powerful 
superconducting magnets are not available.) If one imposes a limit on the maximum depth 
of the decay ring, then there is a maximum baseline for which one can source 3 x 10^*^ ion 
decays per year in the straight section. For all baselines longer than this, the length of 
the straight sections need to be curtailed. If the flux at the Magic Baseline is sufficiently 
reduced, there will become a point when the S cleanliness of the signal is insufficient to 
achieve the best physics reach, even for a single baseline. At such a point, combinations 
involving synergy would then be sought to obtain a broad overall physics reach. By choosing 
to use ^^Ne and ^He only, this study was able to expose both detectors simultaneously to 
the neutrino flux. In studies using all four candidate ions, one effectively cuts the event 
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rate in half since the low-Q and high-Q pairs need to be run separately and irradiate only 
the appropriate detector. The combination of the synergy between baselines and the higher 
event rate from using only two ions is the reason why the relatively minimal configurations 
discussed here have returned similar physics to conservative Magic Baseline proposals. 



6. Summary 

Assessing the physics reach of long baselines with far detectors based at the potential 
LAGUNA sites is of critical importance for strategic decisions towards accelerator and non- 
accelerator neutrino physics in Europe. Recently, as part of the LAGUNA design study, 



work studying the physics return for a Beta Beam with LENA [29| as a far detector, and 



a high powered Superbeam directed towards GLACIER [16| have been performed. In this 
article, this has been continued to include Beta Beam physics with GLACIER, but with 
the mass effectively reduced to 50 kton to offset any shortfall in the useful ion decay rate. 
R&D towards the detector response for the incident beam is currently unavailable, and 
so the approach was to be conservative in energy resolution and efficiency assumptions. 
Configuring the detector with a low energy threshold and binned in a manner to extract the 
oscillatory structure of the signal, the physics for single and double baselines was simulated. 

The physics reach for a Beta Beam directred towards the Slanic Mine [L = 1570 km) 
and Pyhasalmi (L = 2300 km) were simulated initially. The shorter baseline performed 
best owing to its larger event rate and weaker sign(Am|^) degeneracy. Non-zero ^^13 could 
be established down to sin^ 26'i3 ~ 10^'^ for the 7 = 570, 350 facility. The effect of altering 
the boost of the ^^Ne ion is large for the Slanic baseline, especially for CP-violation. This 
indicates that degenerate solutions pose a problem for the less energetic neutrinos. The 
effect was much less apparent for identifying the correct mass ordering. This is because the 
increasing the boost gives access to higher energy neutrinos without substantially increasing 
the flux at low energies. For Pyhasalmi, the sensitivities were much weaker for all indicators 
with little return below sin^ 2^13 = 10~^. Although the matter effect is larger at Pyhasalmi, 
the reduced flux from being more distant from source allows the degenerate solutions to 
remain statistically significant for small 613. 

These results indicated that the data from the Slanic and Pyhasalmi Liquid Argon de- 
tectors should be combined with a concurrent run directed towards a large Water Cerenkov 
based at Canfranc. Such an addition provides an almost clean measurement of 6*13 and 6 
down to sin^ 2^13 ~ 10""^. This information is then available for the longer baseline to re- 
move the degenerate solutions and extract the mass ordering from its data set. The reach on 
the neutrino mass ordering is competitive with the conservative two baseline configurations 
incorporating the Magic Baseline. The access to 613 and CP-violation was dominated by the 
short baseline with reach down to sin^ 2^13 = 5 ■ 10~^. 

In conclusion, the results of this study indicate that European Beta Beams can be pro- 
vide competitive physics reach to move extensive proposals. The optimal Beta Beam con- 
figuration, therefore, has yet to be determined and should be investigated further, taking 
into account available knowledge on technological restrictions. More generally, these results 
demonstrate that the Magic Baseline is not mandatory for determining the sign of Am'^^: 
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the combination of data from a short basehne with data heavy in degeneracy can be an 
equally as powerful phenomenological tool. 
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